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Abstract This paper deals with the study of the reflection and transmission phenomena when axial symmetric 
body waves incident on the base of a poroelastic semi-infinite solid cylinder, surrounded by another medium. 
Cylinder is assumed to be isotropic so that Biot's theory of poroelasticity can be employed. Reflection and 
transmission coefficients are computed as a function of angle of incidence in the case of permeable base. In 
addition, square root of energy ratio is computed for the body waves. Numerical results are presented graphically 
for two types of poroelastic solids, namely sandstone saturated with kerosene and sandstone saturated with 
water. 
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1 Introduction 


The wave reflection and transmission phenomena have long history and received much attention during the 
past half-century in many scientific fields, such as Marine Seismology, Geotechnical Engineering, Acoustics 
and Geophysics. Analysis of reflection and transmission phenomena can be used to understand the various 
materials. Sandstone is a great source for quartz which is very useful in day-to-day life. Sandstone deposits 
are cylindrical in shape and surrounded by poroelastic medium soil or rock. The study of reflection and 
transmission of waves incidented at the interface gives the information pertaining to sandstone deposits. The 
reflection and transmission coefficients are generally based on the medium consisting of two non-homogeneous 
layers separated by a horizontal interface. In Acoustics, the said coefficients are used to understand the effect of 
various materials on their acoustic environments. In the research domains such as Geophysics and Medicine, 
the pertinent analysis can be used as nondestructive evaluation (NDE) tool. Displacement discontinuities 
conserve energy and yield frequency-dependent reflection and transmission coefficients. To the best of author's 
knowledge, some of the contributions in this domain are as follows. Reflection and transmission coefficients 
for harmonic plane waves incident at arbitrary angles upon a plane linear slip interface are computed in the 
case of elastic media [1]. Employing transfer matrix method, Bogy and Gracewski [2] derived the plane wave 
reflection coefficient for a layered solid half-space. In the paper [2], the reflection coefficient is derived for 
an Isotropic, homogeneous elastic layer of arbitrary thickness that is perfectly bonded to such an elastic half- 
space when plane waves are incident from an inviscid fluid onto the layered solid. Singh [3] discussed the 
reflection of P- and SV -waves from the free surface of an elastic solid with generalized thermo-diffusion. 
In the paper [3], it is concluded that the reflection coefficients depend upon the angle of incidence of P- and 
SV-waves and thermo-diffusion parameters. Amplitude ratios and energy ratios are computed numerically, 
and the conservation of energy across the interface is verified for an elastic solid half-space [4]. Kumar et al. 
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[5] investigated the reflection and transmission of plane waves at the loosely bonded interface of an elastic 
solid half-space and a micropolar thermo-elastic diffusion half-space. In the paper [5], the effects of relaxation 
times and loosely bonding parameters shown in energy ratios for a specific model. The analytical expressions 
for amplitude ratios, reflection and transmission coefficients in terms of angles and material parameters have 
been obtained for quasi-longitudinal and quasi-transverse wave incidence [6]. Reflection and transmission of 
elastic waves from the boundary of a fluid-saturated porous solid are discussed by Wu et al. [7]. Santos et 
al. [8] showed the importance of inclusion of the frequency correction factor to analyze wave propagation 
for frequency varying between the seismic and ultrasonic ranges. Propagator matrix method is applied for 
calculation of reflected coefficients in the case of inhomogeneous, anisotropic, poroelastic seafloor [9]. In 
the paper [9], the matrix is integrated using an implicit technique. For generalized thermo-elastic medium, 
the reflection phenomena of SV-waves are discussed by Abd-Alla and Al-Dawy [10]. In the paper [10], the 
reflection coefficient is calculated for shear wave that incident from within the solid on its boundary. Lin et 
al. [11] studied the reflection of plane waves in a poroelastic half-space saturated with inviscid fluid. The said 
paper reports the amount of solid frame stiffness controls the response of a fluid-saturated porous system. 
Employing Biot’s theory of poroelasticity [12], the reflection of plane waves at boundaries of a poroelastic 
half-space is investigated by Tajuddin and Hussaini [13]. In the paper [13], it is clear that the overlapping 
parameter between solid and fluid plays a significant role in generating the reflected slow dilatational wave. By 
taking the sealed pore and the open pore boundaries, wave reflection and transmission are discussed in [14]. 
The said paper revealed that the amplitudes of the reflected and transmitted waves strongly depend on boundary 
conditions. Novel expressions for the reflection and transmission coefficients in terms of the spectral properties 
of the governing differential system are given in [15]. Cui and Wang [16] studied reflection and transmission 
of plane waves at an ocean floor interface. The results of this paper [16] indicate that the effect of the squirt 
flow is noticeable in the reflection and refraction phenomena. Semi-relativistic reflection and transmission 
coefficients for two spinless particles separated by exponential and rectangular-shaped potential barriers are 
investigated by Ethyl et al. [17]. The reflection and transmission of plane waves between two different fluids 
saturated with porous half-spaces are given in [18]. In the paper [18], the variations of amplitude ratios with 
angle of incidence are calculated. Corredor et al. [19] investigated scattering problem at all angles of incidence 
for a single layer embedded between two half-spaces with dissimilar media. Most of the works in this domain 
are confined to half-spaces. However, to the best of author's knowledge, study of reflection and transmission 
coefficients is not made extensively in the case of poroelastic cylindrical solids. The propagation of thick- 
walled hollow poroelastic cylinder under plane strain conditions is studied by Reddy and Tajuddin [20]. The 
said paper describes transitions from the case of plate to thin shell, and to thick-walled hollow cylinder, and then 
to solid cylinder. Investigation of reflection and transmission of acoustic waves of a porous solid is made using 
experimental reflected and transmitted signals [21]. Motivated from the results of said papers, reflection and 
transmission phenomena are analyzed in the case of a poroelastic solid cylinder. In the present work, reflection, 
transmission coefficients and energy ratio are computed for the fast P-wave (P — I), the slow P-wave (P — IT) 
and the shear (SV)-wave in the case of permeable base. 

This paper is organized as follows. In Sect. 2, formulation and solution of the problem are given. Boundary 
conditions are given in Sect. 3. In Sect. 4, numerical results are described. Finally, conclusion is given in 
Sect. 5. 


2 Formulation and solution of the problem 
Consider an isotropic poroelastic solid cylinder whose axis is along z-axis of cylindrical coordinate system, 
surrounded by another medium. Assume axial symmetric body waves that incident at z = 0 plane. Because of 


medium change, two phenomena reflection and transmission do occur. The geometry of the problem is given in 
Fig. 1. The equations of motion of a homogeneous, isotropic poroelastic solid in absence of dissipation are [12]: 


* ü^ , " 
NV?ii + V(Ae + Qe) = ari ule + pU), 
a? " 
V(Qe+ Re) = ag 024 + p20). (1) 
where 4 = (u,v, w) and U = (U, V, W) are the displacement vectors of solid and fluid, respectively, e 


and e are the dilatations of solid and fluid, and P = A + 2N, Q and R are all poroelastic constants. The 
mass coefficients p11, 012 and p22 satisfy the relations pj = p11 + p12 = (1 — B)ps, po = P12 + 222 = Boer, 
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Fig. 1 Geometry of the problem 


p = pi p» = Ps + (pf — ps) [13]. Here, p, and p2 are mass of solid and fluid per unit volume of aggregate, 
p is the total mass of solid-fluid aggregate per unit volume, p11, 012 and p22 are the mass coefficients, o, and 
pr are the mass density of the solid and fluid, and £ is the porosity of the aggregate. The solid stresses o;; and 
fluid pressure s are given by [12] 


oj; = 2Neij + (Ae + Q&e)üjj, (i, j—1,2,3, s= Qe+ Re. (2) 
In Eq. (2), 5;; is the well-known Kronecker delta function. 


Consider the following potential decomposition of displacements in the case of axial symmetric waves 
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The displacement potentials $'s and y’s functions of r, z and t are assumed as given below. 
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In the above, o; is the angle made by the kth incident wave with the normal, a —angle made by the 
Ith reflected wave with the normal when kth wave is incident, $(; —incident wave potential, p% —reflected 
wave potential, o—frequency, 6,—wavenumber of incident wave, 0j —wavenumber of reflected wave, and 


Z M = tam / Ax)—Treflection coefficient of Ith reflected wave when kth wave is incident. Using Eq. (4) in (3), 
after a long calculation, the following solid displacements are obtained. 
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2 
—i[ókó COS oy + 2 5g” cos a + ($49 Jr) + 53g sin o5 (Dg = = 1,2. 
w= (5) 
{debe sin o + > 5g” cos af? + (6° Jr) + 539 sin a, (y = —3. 


Using these displacement components in Eq. (2), the nonzero stress components and fluid pressure are 
derived and are given by 
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In Eq. (6), yy = un E Y, Ak = inu LU E f , k,l = 1,2, 3, where Vi and 


V» are the velocities of dilatational waves of the first and second kind, respectively, and V3 is the velocity of 
shear wave. 


3 Boundary conditions 


For permeable and stress-free base, the boundary conditions are given by 


Ozz +s = 0, orz = 0, s = 0, atz = 0. (7) 
From Eqs. (6) and (7), the following relations are obtained. 
: k) = (k) k 
(k) (k) 
>a; Z“ = Aio k,l = 1,2,3, (8) 
l=1 
where AT are 
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AT = (N + ((A + Q) + (Q  R)uj) sin 2o,)02 —2((A+ Q) + (Q + R)u cos aj )ó3), 
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AT = (Q + Ru) cos2o43? — 2(Q + Rui) sinorój, k = 1,2, m = —(Q + Ru3) cos œ3ô3, 
AQ) = (-2N cos? ax + (A + Q) + (Q + Ruy) cos 24x)? 
—2((A+ Q) + (Q + R)ucsinogóg), k= 1,2, 
A = Nsin 2082, k= 1,2, AQ) — N cos 2063, 
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When a fast P-wave (P — I) is incident, slow P-wave (P — II) is incident and for a shear SV -wave, the 
reflection coefficients are 


mE n Y vk = 1923; 


(k) 
Zi 

(10) 

Similarly, the transmission coefficients are computed for body waves as that of reflection case. The expres- 


sions for transmission coefficients are similar as those of Eqs. (8)-(10), where Cos function is replaced by Sin 
function. 


4 Numerical results 


The reflection, transmission coefficients and energy ratio are computed for permeable base by introducing the 
non-dimensional parameters given by 
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are the normalized phase velocities of kth incidented and /th reflected waves, im In the above, 


Vom ,/ a is the reference velocity, V; and V2 are the velocities of dilatational waves of the first and second 


kind, respectively, and V3 is the velocity of shear wave. First reflection coefficient Z oo is computed and hereby 
energy ratio is computed using Eq. (12) for permeable base. For illustration purpose, two types of poroelastic 
solid cylinders are used, namely cylinder-I made up of sandstone saturated with kerosene [22] and cylinder-II 
made up of sandstone saturated with water [23]. The physical parameter values and normalized phase veloc- 
ities of these materials are given in Tables 1 and 2. Employing these values in Eq. (9), reflection, transmission 
coefficients and energy ratios are computed as a function of angle of incidence for a fixed wavenumber and the 
results are depicted in Figs. 2, 3, 4, 5, 6, 7, 8,9, 10, 11, 12, 13, 14, 15, 16, 17, 18 and 19. The notations RCCyI-I, 
RCOCyI-IL TCCyl-I, TCCyl-II, ERCyl-I and ERCyI-II used in the figures represent the reflection coefficient of 
cylinder-I, reflection coefficient of cylinder-II, transmission coefficient of cylinder-I, transmission coefficient 
of cylinder-II, energy ratio of cylinder-I and energy ratio of cylinder-II, respectively. The following particular 
cases are considered. 
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Table 1 Material parameters 


Material parameters aj a2 a3 a4 dj d» da 
Cylinder-I 0.843 0.065 0.027 0.234 0.901 —0.001 0.1 
Cylinder-II 0.96 0.006 0.029 0.412 0.876 0 0.124 


Table 2 Normalized phase velocities 


1 1 


Material YM a2 VA3 
Cylinder-I 1.0008427 0.4531887 0.5096217 
Cylinder-II 1.047029 0.4832048 0.6857985 
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Fig. 2 Reflection and transmission coefficients as a function of angle of incidence (in degrees) for Z d 
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Fig. 3 Reflection and transmission coefficients as a function of angle of incidence (in degrees) for a 
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4.] Incident fast dilatational wave 


(1) (1) 


The reflection coefficient Z; (i = 1,2, 3) and transmission coefficient T; (i = 1,2,3) are obtained as a 
function of angle of incidence in the case of P — I wave, and the results are depicted in Figs. 2, 3 and 4. From 
these figures, it is seen that as angle of incidence increases, the reflection coefficient, in general, increases 
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Fig. 5 Energy ratios as a function of angle of incidence (in degrees) for E n 


B 5 
o 45 
c 4 
9 35 —e— ERCyH 
Oo 
38 25 —E-— ERCyI-II 
of 2 

1.5 
$ 1 
5 0.5 
e 0 

5 10 15 20 25 30 
Angle of incidence 
Fig. 6 Energy ratios as a function of angle of incidence (in degrees) for E 

> 
S 26 
2 —e— EROyH 
o 21 
5 2 16 —H- EROyHI 
oo 11 
o 6 
S 1 
3 
e 5 10 15 20 25 30 


Angle of incidence 


Fig. 7 Energy ratios as a function of angle of incidence (in degrees) for E 
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Fig. 8 Reflection and transmission coefficients as a function of angle of incidence (in degrees) for Z a 


and transmission coefficient, in general, decreases for both the cylinders. In the case of cylinder-II, reflection 
coefficients Z D, £z zi are negative. When the angle of incidence is 15°, reflection coefficient (Z aa has 


maximum value. When the angle of incidence is 25°, transmission coefficient (T®) has a minimum value. The 
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Fig. 9 Reflection and transmission coefficients as a function of angle of incidence (in degrees) for ze 
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Fig. 10 Reflection and transmission coefficients as a function of angle of incidence (in degrees) for ge 
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Fig. 11 Energy ratios as a function of angle of incidence (in degrees) for E 2) 
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Fig. 12 Energy ratios as a function of angle of incidence (in degrees) for ED) 


energy ratio E a) (i = 1, 2, 3) for fast (P — I) wave, slow (P — II) wave and shear (SV) wave is computed 
as a function of angle of incidence, and the results are depicted in Figs. 5, 6 and 7. In Fig. 5, itis seen that the 
energy ratios for cylinder-I are, in general, lower than that of cylinder-II. From Figs. 6 and 7, it is clear that the 
energy ratios for cylinder-I are, in general, higher than that of cylinder-II. From Fig. 7, it is clear that energy 


ratios (EP ) are very high in the case of reflected shear wave. 
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Fig. 14 Reflection and transmission coefficients as a function of angle of incidence (in degrees) for Z D 
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Fig. 15 Reflection and transmission coefficients as a function of angle of incidence (in degrees) for ZO 
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Fig. 16 Reflection and transmission coefficients as a function of angle of incidence (in degrees) for e 


4.2 Incident slow dilatational wave 


The reflection coefficients 2? — 1, 2, 3) and transmission coefficients Te (i — 1,2, 3) are obtained as a 
function of angle of incidence in the case of slow P — II wave. The pertinent results are depicted in Figs. 8, 9 
and 10. From these figures, it is seen that as angle of incidence increases, the reflection coefficient, in general, 


increases, but the transmission coefficient, in general, decreases. The reflection coefficients Z n. zi ; zo 
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Fig. 19 Energy ratios as a function of angle of incidence (in degrees) for E 


are negative for cylinder-II. When angle of incidence is 25°, the reflection coefficient has maximum value and 
has minimum value when the incident angle is 30?. The energy ratios are computed against angle of incidence, 


and the results are depicted in Figs. 11, 12 and 13. From these figures, it is clear that E e gradually increases. 
Also, it is clear that energy ratios for cylinder-I are, in general, higher than that of the cylinder-II. From Fig. 13, 


it 1s observed that energy ratios (EP) are very high in the case of reflected shear wave. When the angle of 
incidence is 15°, energy ratios (EP) of both the cylinders coincide in the case of reflected shear wave. 


4.3 Shear incident wave 


The reflection coefficients z® (i = 1, 2, 3) and transmission coefficients pe (i — 1, 2, 3) are computed as 
a function of angle of incidence in the case of shear SV wave. The relevant results are depicted in Figs. 14, 
15 and 16. From these figures, it is clear that the reflection coefficient, in general, increases in the case of 
cylinder-I, and the result is reversed in the case of cylinder-II, whereas the transmission coefficient, in general, 
decreases for both the cylinders. When the angle of incidence is 30°, the reflection coefficient of cylinder-II 


has minimum value. Energy ratios E E ; ED ,E 9 are computed for fast (P — I) wave, slow (P — IT) wave 
and shear (SV) wave, and the results are depicted in Figs. 17, 18 and 19. From these figures, it is observed that 
energy ratios of cylinder-II are higher than that of cylinder-I. From Fig. 17, it is observed that energy ratios of 
cylinder-I are slightly linear in nature. 
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5 Conclusion 


Employing the Biot’s poroelasticity theory, two wave phenomena, namely reflection and transmission, are 
investigated. The energy ratios are computed for two dilatational waves and one shear wave in the case of 
permeable base. For numerical process, two types of solids , namely sandstone saturated with kerosene (say 
cylinder-I) and sandstone saturated with water (say cylinder-II), are employed, and the results are presented 
graphically. From the numerical results, it is observed that in the case of cylinder-I, reflection coefficient for fast 
and slow dilatational waves increases and, for shear wave, it decreases. In the case of cylinder-II, the results are 
reversed. Both the materials are sandstone related and differ only in the fluid part. Thus, it can be inferred that the 
fluid part causes the above discrepancy. For all the waves, reflection and transmission coefficients for cylinder-II 
are negative. The energy ratios increase for all the waves, in the case of both cylinder-I and cylinder-II. 
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